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Abstract 

We consider a hypothetical apparatus that implements measure- 
ments for arbitrary 4-local quantum observables A on n qubits. The 
apparatus implements the "measurement algorithm" after receiving 
a classical description of A. We show that a few precise measure- 
ments, applied to a basis state would provide a probabilistic solution 
of PSPACE problems. The error probability decreases exponentially 
with the number of runs if the measurement accuracy is of the order 
of the spectral gaps of A. 

Moreover, every decision problem which can be solved on a quan- 
tum computer in T time steps can be encoded into a 4-local observable 
such that the solution requires only measurements of accuracy 0(1/T). 

Provided that BQPt^PSPACE, our result shows that efficient algo- 
rithms for precise measurements of general 4-local observables cannot 
exist. We conjecture that the class of physically existing interactions 
is large enough to allow the conclusion that precise energy measure- 
ments for general many-particle systems require control algorithms 
with high complexity. 



* e-mail: {woe j an , j anzing , decker}@ira . uka . de 
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1 Measuring /c-local observables 



A characteristic feature of quantum theory is that there exists an 
abundance of mutually incompatible observables (described by self- 
adjoint operators A) for every quantum system and it is by no means 
obvioTis how to implement measurement procedures for all these quan- 
tities. On a quantum computer one could in principle measure every 
observable A as follows: Find a unitary transformation U which diago- 
nalizes A with respect to the computational basis. Then A is measured 
by implementing U and measuring the logical state of each qubit. By 
identifying each binary word with the corresponding eigenvalue of A 
this procedure reproduces all probabilities 

Pj ■= tripPj) 

correctly where p is the density matrix of the quantum register and 
(Pj) is the family of spectral projections of A. 

However, the implementation of the diagonalizing operation U will 
in general be hard. Therefore, one may restrict the attention to specific 
classes of observables. It is natural to consider observables with physi- 
cal relevance. For example, the qiiantum observable "energy" , mathe- 
matically described by the self-adjoint operator H (the Hamiltonian) , 
is certainly one of the most important observables in physics. It deter- 
mines the dynamical and thermodynamic behavior of the considered 
quantum system. Furthermore the eigenstates of the Hamiltonian, the 
energy levels, are "directly" observable in many physical situations. 
For instance, in spectroscopy the eigenvalues of the Hamiltonian deter- 
mine the frequencies of emitted or absorbed photons. Nevertheless, 
the determination of the energy levels in interacting many-particle 
systems is in general a difficult task. 

To explain this more explicitly, we need to describe the class of 
operators which is considered. First we note that physical interaction 
Hamiltonians usually satisfy some locality condition in the following 
sense. We call an n-qubit operator fc-local if it is a sum of opera- 
tors which act on at most A; particles non-trivially. For fundamen- 
tal interactions between real physical particles one has more specific 
statements and may restrict the attention to pair- interactions. Nev- 
ertheless, fc-local interactions among qubits are physically reasonable. 
They may describe effective Hamiltonians and there is not necessar- 
ily a one-to-one correspondence between qubits and physical particles 
{I qubits may, for instance, describe the state of one particle). The 



2 



following results show that it is in general difficult to compute the 
spectrum of fc-local Hamiltonians. 

The problem of determining the lowest energy value of a (classical) 
spin-spin interaction of Ising type is known to be NP-complete ^Hj- 
For interacting qubits determing the lowest energy value is even QMA- 
complete ( "Quantum-NP" ) if one allows 3-local interactions only 
Note that in these NP and Quantum-NP problems the task is not to 
determine the lowest eigenvalues with high precision. The demanded 
accuracy is only inverse polynomially in the number n of interacting 
qubits. This has implications for the measurement procedure above: 

The unitary U that maps the eigenvectors of A to the computa- 
tional basis states is only helpful for measuring A if the correspondence 
between computational basis states and the eigenvalues of A is known. 
Therefore, this method would require to know the spectrum of A. For 
2-local or 3-local observables one would need the solution of NP- and 
QMA-hard problems, respectively. 

In this paragraph we will explain that measurements of /c-local 
observables A are possible up to inverse polynomial accuracy without 
using any knowledge on the spectrum of A. 

Here we do not need a precise definition of accuracy, we only de- 
mand that the following condition is satisfied: 

Postulate 1 (Measurement accuracy) 

A measurement with accuracy AX has the following property: For all 
density matrices p the probability to obtain an outcome in the interval 
I := [Xj - AX, Xj + AX] is at least (3/4) tv{pPj). 

Our result is not sensitive to the particular definition of accuracy. 
However, it is convenient to work with the formulation above. 

Now we describe how to implement approximative measurements. 
The idea is that for every /c-local A (with k constant) the correspond- 
ing time evolution Ut := exp{—iAt) (if A is interpreted as a Hamil- 
tonian H oi a quantum system) can be simulated efficiently in an 
approximative sense. Explicitly, it has been shown that the simula- 
tion of Ut with elementary gates up to an error of e (with respect to 
the operator norm) requires 0(t^/e) gates [3]. Now we can choose t 
in such a way that there is a one-to-one correspondence between the 
eigenvalues of Ut and A. This is the case whenever < vr. An up- 
per bound on the norm of a /c-local operator is easy to get. We assume 
without loss of generality that each A;-local term is upper bounded by 
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Figure 1: Circuit for performing measurements of an observable A. The 
powers of [/ := exp{—iAt) are implemented as conditional gates controlled 
by the ancilla register. 



the value 1. There are at most 

fc-local terms. Therefore, one has ||^|| = O(n^). For an appropriate 
value of t we can implement measurements of the "observable" Ut 
using^ the quantum state estimation procedure [HI- We will briefly 
sketch the idea. In the following we drop the index t. 

The circuit for phase estimation is shown in Fig. ^ It acts on the 
registers Rs and Ra- Controlled-f/-' gates are implemented in such a 
way that is performed on Rs if and only if the ancilla register Ra 
is in a state corresponding to the binary word j. This can be done by 
implementing controlled 

gates which are applied if and only if the /th qubit of the ancilla register 
is in the state The algorithm starts with an equally weighted 
superposition 

i=l 

^Notc that also a unitary operator defines in a canonical way an observable by its 
spectral projections if one allows complex measurement outcomes. 
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of all ancilla register states where m is the number of ancilla qubits. 
After applying the controlled [/■'-operation the discrete Fourier trans- 
form of size 2™ is applied to the ancilla register. Then the eigen- 
values of U can be readout with an error of the order 1/2™. Now 
we consider the running time of this scheme (depending on the accu- 
racy). Obviously, this depends on the running time for implement- 
ing the controlled- f/-' operations. The natural method to implement 

= exp(—iAjt) is to simulate the time evolution with Hamiltonian 
A for the time tj. The substitution of the corresponding gates by con- 
trolled gates is straightforward. However, j grows exponentially with 
m. Consequently, this method requires exponential running time for 
exponential accuracy. It is likely that all possible schemes for measur- 
ing A precisely share this disadvantage. 

For a black-box unitary U it is clear that exponential accuracy 
requires exponential time since the black-box unitary U has to be ap- 
plied an exponential number of times. This can be seen by subjecting 
two state vectors to different unitaries U and U with the same eigen- 
states but slightly different eigenvalues. Then [/' and W can only 
lead to distinguishable states for large /. In it is described how to 
convert the time evolution exp{—iHt) according to an unknown pair- 
interaction Hamiltonian if to a controlled-exp(— ii/t) evolution. This 
shows that black-box settings for unknown A do in principle make 
sense for energy measurements. 

For non-black box interactions we cannot obtain lower bounds on 
the measurement complexity by similar arguments since the appara- 
tus receives a classical description of the observable to be measured. 
However, the result of this paper suggests that even if the interaction 
is known there is no efficient measurement scheme with exponential 
accuracy. We show that measurements of 4-local n-qubit observables 
A could be used to solve PSPACE-problems in polynomial time pro- 
vided that the accuracy is sufficient to distinguish between the differ- 
ent eigenvalues of A. 

One may ask whether there may be any physical processes for 
measuring 4-local observables that do not rely on quantum circuits 
(consisting of elementary gates). For instance, one may guess that a 
measurement of the energy of a system is simpler than a measurement 
of an arbitrary /c-local observable because energy is a specific observ- 
able determining many physical aspects of the system. But note that 
the quantum version of the Strong Church- Turing Thesis (compare 
[HI E] ) states that every problem that can be solved efficiently using 
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some physical process can be solved efficiently by a quantum com- 
puter. 

If there existed any efficient scheme for precise measurements of 
4-local observables this would imply either of the following statements: 

1. The measurement process cannot efficiently be simulated on a 
quantum computer (in contrast to the strong quantum Church- 
Turing Thesis). 

2. There are polynomial time algorithms to solve probabilistically 
PSPACE problems, i.e., PSPACE = BQP. 

Assuming that both implications are unlikely, our result strongly sug- 
gests limitations for future quantum measurement technology. 

The structure of the paper is as follows. In the next section we 
consider a class of quantum circuits with polynomial size. They could 
solve PSPACE problems if they were applied an exponential number 
of times on a polynomial number m of qubits. 

In Section 3 we describe how to construct a 4-local observable 
corresponding to this circuit in such a way that precise measurements 
would solve PSPACE problems in polynomial time. 

2 Characterizing PSPACE by circuits 

The complexity class PSPACE is usually defined with respect to the 
Turing machine model PSPACE is the class of all languages rec- 
ognizable by polynomial space bounded deterministic Turing machines 
that halt on all inputs [111. 

For our purposes we need a characterization of PSPACE with re- 
spect to quantum circuits. In particular, we need the result that every 
PSPACE language can be recognized by applying an appropriate cir- 
cuit many times. 

Theorem 1 (PSPACE) 

For every language L in PSPACE there is a polynomial-time uniformly 
generated family of quantum circuits {Vi)i£n consisting of Si = poly{l) 
elementary quantum gates and acting on mi = poly{l) many qubits. 
The circuit Vi decides whether an input string x of length I is an 
element of L in the following sense. 
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There is a polynomial-time computable natural number ri such that 
the ri-fold concatenation ofVi solves the corresponding PSPACE prob- 
lem, i.e. 

Vp {\x) \y) |00 . . . 0)) = \x) |y © f{x)) |00 . . . 0) , 

where f is the characteristic function of L, i.e., f{x) = 1 if x £ L 
and f{x) = otherwise. The vector \x) is the basis state given by the 
binary word x E {0, 1}', the vector \y) is the state of the output qubit 
and |00 . . . 0) is the initial state of mi — I — 1 ancilla qubits. 

Proof. In order to construct the circuit Vi corresponding to a PSPACE 
problem we need to have an upper bound for the required space. This 
is, for instance, the case for the PSPACE-complete problem QBF 
(Quantified-Boolean Formulas). It can be solved within the space 
0(/^) where I is the length of the input. This space bound determines 
mi, the number of qubits. 

Let M be a Turing machine that solves QBF within space 0{l^). 
Now we construct a quantum circuit Vi that simulates the Turing 
machine M for input length /. Since the computational steps of a 
quantum circuit are unitary (thus reversible), we have to work with 
a reversible Turing machine R instead of M (the latter could be irre- 
versible). Each application of the constructed circuit simulates one or 
two steps of R. 

Due to a result of Lange et al. (Theorem 3.3 in J^) it is possible 
to simulate irreversible Turing machines by reversible ones without in- 
creasing the necessary space too much. More precisely, they give the 
simulation of a space-bounded Turing machine M by a reversible Tur- 
ing machine R operating on the same space. In general, the reversible 
simulation by R may have an exponential time overhead. The run- 
ning time overhead is not relevant here because we can derive an upper 
bound on the running time of the reversible machine from the required 
number of qubits. In the following we work with the reversible Turing 
machine R. 

The fact that every Turing machine can be simulated efficiently 
by circuits is standard Here we need an explicit construction 

converting the reversible Turing machine into a circuit consisting of 
reversible gates. 

The circuit acts on the following registers: 

1. The register head encodes the internal state of the Turing ma- 
chine. 
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2. The register tape_index stores the current location of the head. 

3. The register ACC is the accumulator (temporary storage). 

4. The register tape corresponds to a sufficiently large region of 
the tape that is required for computation. It consists of cell 1 to 
cell A'^; where A'^; is the space bound corresponding to the input 
length /. 

Each step of the reversible Turing machine of Lange et al. is either 
a moving or read-and- write transition'^. A moving transition has the 
form p — > (g, ±1). That means that in state p the machine makes one 
step to the right (+1) (respectively to the left (—1)) and changes into 
state q without reading or writing any tape cell. A read-and-write 
transition has the form (p, a) — > (g, b) meaning that in the state p the 
machine overwrites the symbol a with the symbol b and changes into 
state q without moving the head. 

Furthermore, in our construction it is determined by the state of 
the head whether the system performs a moving or a read-and-write 
operation (and not by the state of the tape). In other words, the state 
set Q of -R is the disjoint union of a set Q of read-and-write states, a 
set of right-moving states and a set Q'^ of left-moving states. 

Due to reversibility of R the moving transitions can be imple- 
mented as a unitary transformation on the registers head and tape_index. 
A right-moving transition {p G Q^) translates as follows: 

|p)head ® |^)tape_index * |9)head |i + l)tape_index • (1) 

Analogously, a left-moving transition {p G Q'^) translates as follows: 

|p)head ® |^)tape_index * |9)head ® \i — l)tape_index • (2) 

Note that the operations on the register tape_index are computed 
modulo N , where N is the number of tape cells. Although the Turing 
machine will never move to the right when the head is at position 
and never to the left when it is at position 1, this definition guarantees 
that eqs. and ((21) define unitary operators. 

These transformations can be realized efficiently as a unitary trans- 
formation C/moving acting only on the registers head and tape_index. 

Again, due to reversibility of R the read-and-write transitions can 
be realized by a unitary transformation. Let (p, a) — > (g, b) be a read- 
and-write transition. There is a unitary transformation W^jy^ acting 

^This separation is useful in order to characterize reversibility of Turing machines |14j . 
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on the registers head and ACC reaUzing 

|p)head ® |a)ACC —>■ head ^ |*^/ACC • 

We denote by SWAP(ACC,tape [i] ) the unitary operation that swaps 
ACC and the ith ceh of tape. We denote by A,(C/) the controlled 
operation that performs U if and only if tape_index has the value i. 
Now we define Uj./„ as the concatenation 

N N 

Yl Ai (SWAP(ACC, tape [i] )) W,/„ JJ (SWAP(ACC, tape [i] )) . 

i=l i=l 

The transformation U = C/r/w ^moving is the transformation that corre- 
sponds to a moving and/or read-and-write transition of the reversible 
Turing machine R. (Note that if a read-write translation is followed 
by a moving transition then this U performs both transitions.) 

The constructed circuit U does not satisfy all requirements of the 
theorem. The first problem is that we do not know the running time of 
the reversible Turing machine R. Consequently, we do not know how 
many times we have to perform the elementary circuit U to obtain 
the solution. Even if we knew how many times we have to apply U, 
the corresponding transformation would in general change the state 
\x) and produce some garbage on the ancillas. 

To circumvent the first problem we introduce some idle cycles to 
guarantee that the running time is an efficiently computable function 
ri of the input length I. The second problem is solved by uncomputing 
the operations carried out during the computational steps and the idle 
cycles. 

In the step r;/2 (proper computational steps and the idle cycles) 
the solution is copied to the register solution. In the following r;/2 
steps we uncompute the idle cycles and the computational steps. The 
computational steps are uncomputed by applying [/^ corresponding 
to running R backwards. 

Now we construct the quantum circuit V that circumvents both 
problems as explained. Note that we drop the index I in the following. 
The circuit operates on the registers head, ACC, tape_index,tape 
and the new registers solution, operation_mode, idle_counter, and 
counter (See Fig. [21). 

The register operationjnode indicates whether the current oper- 
ation is U, idle cycle, reverse idle cycle, or U'^ . These 4 subroutines of 
the whole circuit can be seen in Fig. |21 
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Figure 2: Quantum circuit satisfying the requirements of Theorem 1. The 
gates INC and DEC increment and decrement the register counter and idle 
counter, respectively. The gates 6i62/&'i&2 swaps the state \b1b2) and |&'i&2)- 
The gate 00/01 10/11 is controlled by g/ and the state of the idle counter. 
The symbol qj represents all final states of the Turing machine R. The 
symbol af denote the solution f{x). The bit-fiip on the register solution is 
controlled by and the state of the register counter. 
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Figure 3: The 4 subroutines of the circuit V. Simultaneously with the sub- 
routines the counter runs forward or backward. 

The content of the register counter is incremented after each appli- 
cation of U or idle cycle and decremented after each application of 
or reverse idle cycle. Our construction uses the following upper bound 
on the number of necessary applications of U. Since the number of ba- 
sis states of the register that U acts on is 2™ it does not make sense to 
have r > 2™. Therefore the counter is incremented until the register 
has the resuh 2*"+^ - 1 in order to ensure that the number of required 
applications of U is exceeded. As soon as this number is reached the 
result of the computation is copied to the register solution, i.e., the 
register is incremented by 1 if and only if the answer is "true" . Prom 
this moment on the counter and the idle_counter are decremented. 
As soon as the idle_counter reaches the operation mode is changed 
such that the concatenated application of W is started. After the ap- 
propriate number of applications the initial state of all registers are 
restored except from the register solution which is incremented by 1 
if and only if the answer of the PSPACE problem is "true" . 

Explicitly, one has the following rules: 

1. operation mode 00: perform the circuit U, increment counter 

2. change operation mode 00 01 if idle_counter is 00 • • • and 
head is in a final state 

3. operation mode 01: increment counter and idle_counter 

4. increment solution if operationjnode is 01, counter is 11 • • • 1 

and the first tape cell is in a state indicating if the answer is true 
(we assume that this tape cell contains the result f{x)) 

5. change operation mode 01 — > 10 if counter is in 11 • • • 1 

6. operation mode 10: decrement counter and idle_counter 
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7. change operation mode 10 11 if idle_counter is in 00 - • - 
and head is in a final state 

8. operation mode 11: perform the circuit U\ decrement counter 
and idle_counter 

9. change operation mode 11 00 if counter is in 00 • • • 

Note that the circuit V has the following property: applied to the 
initial state \x) |0 . . . 0) the orbit length is r = 2(2"*+^ — 1) if the answer 
is "false" and 2r whenever the answer is "true". □ 

The dependence of the orbit length on the solution is essential in 
the following section. 

3 Constructing the observable 

In this section we construct a family of observables {Ai) in such a 
way that the spectral properties of Ai reflect the length of the orbit 
(Vf'lx) |00 . . . 0))jgN for inputs of length I. The idea to construct 
Hamiltonians corresponding to quantum circuits already appeared in 
[15j . In this article, the purpose was to show that a closed quantum 
(Hamiltonian) system can in principle implement a circuit without 
any external control operations. Similar constructions were also used 
in the context of complexity theory in order to show that determining 
the spectrum of physical Hamiltonians may be computationally hard 
[21 However, their constructions deal with quantum circuits of 
polynomial size. The whole sequence of gates is in some sense encoded 
into the Hamiltonian. The solution of a NP or QMA problem is then 
reflected in the least eigenvalue of the Hamiltonian. The fact that the 
determination of the least eigenvalue encompasses NP or QMA even if 
only inverse polynomial accuracy is required is due to the polynomial 
length of the program. Here we have typically an exponential number 
of applications and the solution of the problem is therefore encoded 
in the "hyperfine structure" of the spectrum. 

Let y be a quantum circuit as in Theorem ^ and s be its size, 
i.e., the number of elementary two-qubit gates. We need a register 
clock indicating which gate is applied. It consists of s/ qubits. The 
allowed states of the register clock are of the form |0 • • • 010 • • • 0) 
indicating which gate of V is applied currently. We denote by Vj the 
elementary gates of V (in contrast to the preceeding section where the 
index denoted the input length). 
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We first define the forward-time operator 

F = Fi® |l)2(0|2®|0)i(l|i + 

F2(X)|1)3(0|3«>|0)2(1|2 + 

|l)i(O|i0|O),(l|,. 

The operators Vj operate on all registers of the preceding section. The 
operators |0)i(l|i and |l)i(0|i are annihilation and creation operators, 
respectively, on the ith qubit of the clock. 
We denote the linear span of the vectors 

F^>o) forjGN 

with |*o) := \x) (g) |00 • • • 0) (g) |100 • • • 0) as O. All states of this orbit 
are orthogonal until one has recurrence to the initial state |^'o)- This 
can be seen as follows: If the register clock is in an allowed state 
there is only one summand of F that is relevant. Its action on the 
clock is simple since it moves the 1 to the next qubit. Therefore it is 
clear that the first s — 1 states are orthogonal. The whole circuit V 
is a classical logical operation which permutes basis states. Therefore 
the state F*|^o) is either orthogonal to |^'o) or both states coincide. 
Along the same line we argue that all states of the orbit are orthogonal 
until a state coincides with the initial state. Hence F acts as a cyclic 
shift on O. 

The dimension of O is 2sr if f{x) = 1 and sr if f{x) = 0. We 
denote the dimension by d. 

Let w be a primitive complex d-th root of unity. The eigenvalues 
of F restricted to O are 

UJ ,UJ ,UJ2i ■ ■ ■ 1^ 

Furthermore, the initial state vector |^'o) is a superposition of all 
eigenvectors of F restricted to O with equal weights. All this follows 
from properties of the cyclic shift operator. 

The backward-time operator is defined as the adjoint of F. The 
observable A is defined as the sum of the forward and backward time 
operators, i.e., A := {F + F^)/2. It is 4-local since each Vj is 2-local 
and is coupled to a 2-local propagator. The dynamics of the clock 
may be interpreted as a propagation of a spin-wave. Note that the 
idea of our construction is not to implement the quantum circuit V 
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by the autonomous time evolution exp(— i^t). The aim is rather to 
obtain an observable such that its spectral properties correspond to 
the orbit length of the circuit. 

Since F and commute on O the eigenvalues of A restricted to 

are {uj^ + a)^)/2 = cos{2'k j / d) . The non-real eigenvalues are 2-fold 
degenerated. Only the eigenvalues 1 and —1 have multiplicity 1. In 
a hypothetical energy measurement applied to the initial state vector 
l^o) one would obtain all 2-fold degenerated eigenvalues with proba- 
bility 2/d each and the non-degenerated eigenvalues with probability 
1/d. Note that only the first case is relevant for large d since there 
are at most two non-degenerated values. 

Note that d depends on the solution of the PSPACE problem. 
Explicitly, the possible measurement results are 

1. either 

cos(27rj/(2s r)) , j = 0, . . . , 2s r - 1 

2. or 

cos(27rj/(sr)) , j = 0, ...,sr-l 

depending on whether f{x) = 1 or f{x) = 0. 

Note that a perfect energy measurement can distinguish between 
the two cases even after few samples: after applying the function 
"arccos" we obtain values with distance 2Tr/d and all values occur with 
equal probability (if the non-degenerated values are neglected). Then 
it is easy to distinguish between the two cases d = r s and d = 2r s. 

Now we examine what accuracy is sufficient to distinguish be- 
tween the two cases. For doing so, we will restrict our attention to 
those measurement values which are between 1/V2 and — 1/\/2- This 
means that half of the measurement outcomes have to be ignored be- 
cause the probability to obtain an outcome in this interval is about 
1/2. These values correspond to angles in the interval [it/ 4, Sir/ 4] and 
[57r/4,77r/4]. 

In the following we assume that we have obtained a measurement 
value in this interval. For each outcome E we chose j such that 

1 arccos (£') — 2Trj/{2rs)\ is minimal. If f{x) = 1 then the proba- 
bility of obtaining an odd value for j is at least 1/2 • 3/4 = 3/8. If 
f{x) = then the probability of obtaining an even value for j is at 
least 1 • 3/4 = 3/4. Therefore, the probability of odd value is at most 
1/4. This difference in probability allows to distinguish between the 
two cases. It is obvious that the error probability decreases exponen- 
tially with the number of measurements. 
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Note that the observable A has spectral gaps that are consider- 
ably smaller than the required accuracy. This can already be seen if 
we consider the ^-invariant subspace O. The distance of the largest 
eigenvalue 1 and the second largest eigenvalue cos(27r/(i) of H is ap- 
proximatively given by (27r/d)^ since the derivative of the cosinus func- 
tion at is 0. 

Note that the required accuracy is directly connected with an up- 
per bound on the running time T. In our setting the running time 
is the number r of necessary applications of the circuit V times the 
number s of gates of V. 

In the construction of the preceding section we obtained the upper 
bound on rs from the required space. More generally, whenever we 
know that r applications of V are sufficient we need a measurement 
with accuracy of the order l/(rs) to determine the solution of the 
PSPACE problem. 

This discussion proves the following theorem: 

Theorem 2 (Measurement precision vs. running time) 

Let {Ai} be a family of A-local observables corresponding to the family 
{Vi} of quantum circuits in Theorem 1. Then every measurement in 
the sense of Definition^ could be used to solve PSPACE problems in 
polynomial time whenever the accuracy is of the order of the spectral 
gaps of H. It is even sufficient to have an measurement error 1/Ti, 
where Ti is the running time of the algorithm based on the circuit Vi. 

4 Conclusions 

We have shown that every apparatus which implements precise mea- 
surements of 4-local n-qubit observables would solve PSPACE prob- 
lems. This conclusion does only hold for exponentially small errors 
of the measurement. On the other hand, we have argued that al- 
gorithms which measures with inverse-polynomial accuracy can be 
implemented efficiently. Provided that PSPACE problems cannot be 
solved efficiently, i.e., PSPACE ^ BQP, the complexity of measure- 
ments depend on the required accuracy. The statement that exponen- 
tial accuracy has stronger computational power is also well-known in 
classical analog computational models [THl ITT]. 

One may ask why one should try to measure general 4-local ob- 
servables. A possible motivation to develop a complexity theory of 
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measurements is that some proposals for quantum algorithms use joint 
observables on the quantum register |T^ . 

Another motivation is that one is interested in measurements for 
physically relevant joint observables like energy. The reader may ob- 
ject that the specific interactions constructed in this paper are rather 
unphysical for several reasons: 

1. Most interactions in nature are pair- interactions and not 4- local. 

2. Our construction uses long-range interactions among distant qubit 
quadruples. 

3. Interactions in natural many-particle systems have typically high 
symmetry. For instance, the interactions in solid states respect 
the translational invariance of the lattice. 

4. There exist only a few fundamental interactions in physics. 

We have already argued that pair-interactions between particles 
may correspond to /c-local terms if some qubits encode the physical 
state of one particle. This refutes the first objection. 

We conjecture that the solution of PSPACE problems would even 
be possible if the class of observables was restricted to those which 
appear as Hamiltonians of real many-particle systems. This conjecture 
is supported by the following ideas: 

Quantum cellular automatons like the Hamiltonian dynamical sys- 
tem constructed in |19j are also computationally universal. Hamilto- 
nians for those types of cellular automata have the property that every 
cell interacts only with some cells in its neighborhood. By merging 
some cells together to one cell we can always obtain a Hamiltonian 
with pair- interactions among qudits. This seems to indicate that nei- 
ther the symmetry nor locality assumptions on the interactions pre- 
vents the Hamiltonian from corresponding to computationally univer- 
sal networks. Due to the fact that computers exist it is clear that the 
structure of the fundamental interactions is general enough to allow 
universal systems. Therefore we guess that spectral properties of more 
realistic Hamiltonians encode PSPACE problems in a similar way as 
in our paper. 
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